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Abstract 

We use the continuum random phase approximation to describe the muon 

capture on 12 C, 16 and 40 Ca. We reproduce the experimental total capture 

rates on these nuclei to better than 10% using the free nucleon weak form 

factors and two different residual interactions. However, the calculated rates 

for the same residual interactions are significantly lower than the data if the 

in-medium quenching of the axial- vector coupling constant is employed. 
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The capture of a negative muon from the atomic Is orbit, 

H~ + (Z, N)^u fi +(Z-1,N+ 1)* (1) 

is a semileptonic weak process which has been studied for a long time (see, e.g., the reviews 
by Walecka [[1J or Mukhopadhyay |2j and the earlier references therein). The total capture 
rate has been measured for many nuclei; in some cases the partial capture rates to specific 
states in the daughter nucleus have been determined as well. 

The nuclear response in muon capture is governed by the momentum transfer which is 
of the order of muon mass. The energy transferred to the nucleus is restricted from below 
by the mass difference of the initial and final nuclei, and from above by the muon mass. 
The phase space and the nuclear response favor lower nuclear excitation energies, thus the 
nuclear states in the giant resonance region dominate. 

Since the experimental data are quite accurate, and the theoretical techniques of eval- 
uating the nuclear response in the relevant regime are well developed, it is worthwhile to 
see to what extent the capture rates are understood and, based on such comparison, what 
can one say about possible in-medium renormalization of the various coupling constants. In 
particular, there are various indications that the axial- vector coupling constant in nuclear 
medium is reduced from its free nucleon value of qa = 1-25 to the value of qa — 1- The 
evidence for such a renormalization comes primarily from the analysis of beta decay between 
low-lying states of the (sd) shell nuclei . In addition, the "missing Gamow- Teller strength" 
problem, as revealed in the interpretation of the forward angle (p, n) charge-exchange reac- 
tions m, is also often quoted as evidence for quenching of qa- The Gamow- Teller strength is 
concentrated in the giant GT resonance at excitation energies not very far from the energies 
involved in the muon capture, although this latter process is dominated by the transitions 
to the negative parity spin-dipole states. 

Muon capture also depends on the induced pseudoscalar hadronic weak current. At the 
free nucleon level the corresponding coupling constant is determined by the Goldberger- 
Treiman relation [g| 
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= 2M^ 

where m n is the pion mass and Fa(0) = qa = 1.25. (In muon capture one often uses a 
dimensionless quantity gp = m^Fp^q 2 ) at the relevant momentum transfer q 2 ~ —0.9m 2 , 
such that gp ~ 8.4 for free protons.) In nuclear medium Fp can be again renormalized, and 
this renormalization does not necessarily obey the Goldberger-Treiman relation ||. 

The Continuum Random Phase Approximation (RPA) has been used successfully in the 
description of the nuclear response to weak and electromagnetic probes @. The method 
combines the usual RPA treatment with the correct description of the continuum nucleon 
decay channel. We have used this method for calculations of the muon capture processes on 
12 C, 16 0, and 40 Ca. As residual interactions we adopted the finite-range force |§ based on 
the Bonn potential || and the zero-range Landau-Migdal force with the parametrizations 
for 12 C and 16 taken from Refs. |H| and [11], respectively. For 40 Ca we used the standard 



parametrization as, for example, given in Ref . [12] . Note that none of these forces have been 
adjusted to weak interaction data for these nuclei. In the calculation we evaluate the capture 
rate at each energy transfer uj for each multipole separately. The momentum conservation 
condition is fulfilled throughout. We also use an accurate relativistic description of the initial 
bound muon. Form factors and their g-dependence have been adopted from Ref. [[0| . 

The results of our calculations are summarized in Figs. 1-3, which show the capture rate 
as a function of excitation energy, and in Table I, which compares the total capture rates 
with data. The total muon capture rates for 12 C and 16 O, as given in Table I, are defined 
as the part of the rate where the nucleus in the final state is excited above the particle 
threshold and therefore decays via particle emission. The experimental entries in Table I 
were derived by subtracting the partial muon capture rates into the particle bound states 
(we used the average of the various data sets) from the total capture rates. We would like to 
point out that most of the capture rate goes to particle-unbound states. In 12 C capture to 
particle-bound levels (mainly the 12 B ground state) contributes about 16% of the total rate, 
while in 16 the bound-state contributions are roughly 10%. As is obvious from Table I, 



our calculations reproduce the total muon capture rates into the continuum states very well. 
For all three nuclei the Bonn potential slightly overestimates the data, however, by less than 
10%. The Landau-Migdal force reproduces the continuum data for 12 C and 16 remarkably 
well, while it underestimates the 40 Ca data by about 9%. 

For 12 C and 16 partial muon capture rates to particular bound levels have been also 
measured. We compare these data with our calculation in Table II. The calculated partial 
muon capture rates, as given in Table II, again reproduce the magnitude of these rates well, 
with the notable exception of the transition to the 12 B ground state. It is well known that a 
proper description of this Gamow- Teller transition requires additional configuration mixing 
within the p-shell other than provided by the (Ip-lh) RPA approach jrj|]. For the capture to 



the 0~ and 1~ states in 16 N we performed the calculation not only by the continuum RPA, 
but also by the standard RPA which treats all states as bound; the two methods agree with 
each other quite well. For the 2~ state we used only the latter method since the continuum 
RPA gives a much too narrow resonance in this case and the round-off errors are too severe. 

Figure 1 shows the total 12 C(/i~, z^) 12 B capture rate as a function of neutron energy 
above the n + n B threshold in 12 B. Most of this rate goes via 1~ and 2~ multipole 
excitations to the giant dipole and spin-dipole resonances. This is demonstrated in Fig. 1, 
which, as additional information, shows the partial contributions of these multipoles. Giant 
dipole and spin-dipole excitations also dominate the capture rates for the other two nuclei, 
16 and 40 Ca. For example, we find that the \~ and 2~ multipoles together contribute 
about 75% [65%] to the total 16 0(/i", z/ M ) 16 N [ 40 Ca(/i~, ^) 40 K] rate. Excitation of the giant 
quadrupole resonance at about 20 MeV (see Fig. 1) contributes a few percent. For 40 Ca, 
the 0~ and 3~ multipoles each contribute about 10% to the rate. 

In accordance with our discussion above, we find an average excitation energy which in 
all cases corresponds to the regime of the giant dipole and spin-dipole resonances in these 
nuclei. If we consider the Q-values of the reactions, these average energies indicate that the 
average neutrino energy after the capture process is (E v ) ~ 80 MeV, while the remaining 
25 MeV of the muon mass are transferred, on average, to internal nuclear degrees of freedom. 



In Fig. 2, we compare the excitation spectrum for the 16 0(/i~, z/ M ) 16 N reaction as calcu- 
lated for the two residual interactions we used. Both spectra are very similar, but the Bonn 
potential predicts a slightly higher excitation rate to the giant dipole resonances, which ac- 
counts for the 10% difference in the total rates for the two interactions. Finally, in Fig. 3 
we show the capture rate for 40 Ca(/i~, z/ /J ) 40 K as a function of the excitation energy in the 
final nucleus 40 K, i.e., for the whole range of the nuclear response. 

All results presented so far have been obtained using the free nucleon form factors. To 
test the dependence of the calculated rates we repeated our muon capture calculations, 
however, using renormalized values for g A and gp. We present the results in Table III for 
three models of the in-medium renormalization. In model 1 we simply quench the axial 
vector coupling, and keep the Goldberger-Treiman relation intact: 

2m il M p g A 

9 P = — 2 2 ' 9a = 1-0 . (3) 

In models 2 and 3 we modify the relation between gp and g A , using the relation g A = 
fn9nNN I 'M p and two alternative prescriptions given in Ref. |J. So in model 2 

^ m fj,fir9nNN ~ , 1 n /i ok ~ 1 n t a\ 

9 P = : o K / o 2^ n Q^-2 ' 9nNN 9nNN = 1.0/1.25 , g A = 1.0 , (4) 



and in model 3 



. Ql7l n J-kQ-k nn 

g p = —J 5 n -7-2 ' 9-kNN 9-kNN = 1.0/1.25 , g A = 1.0 . (5) 

— q z — 0.7q 1 



As expected from the axial vector dominance of the muon capture cross section, the rates 
in all three cases displayed in Table III are significantly smaller than the ones obtained for 
the free nucleon form factors. Depending on the adopted model for the renormalization, the 
rates are lowered by 20-30%, where the decrease is nearly the same for all three nuclei for a 
given model. This decrease clearly shows that the muon capture rate, although dominated 
by the axial vector current interaction, is also affected by the vector-axial vector interference, 
and by terms containing g P . The induced pseudoscalar coupling decreases the capture rate, 
as one can see best in model 3 where gp is strongly reduced. Comparing the three models 



of renormalization we confirm the known weak sensitivity of the total capture rate to the 
variations in gp. However, most importantly, the rates obtained with the quenched form 
factors are in disagreement with the data. 

To test the sensitivity of our results to the adopted residual interaction we have performed 
calculations for all three nuclei, in which the overall strength of the Landau-Migdal force 
has been decreased or increased by 10%. In all of the calculations qa was set to qa = 1.25. 
Noting that the residual interaction is repulsive in the isovector channel and pushes the 
T = 1 strength to higher excitation energies, a weakening of the force results in the T = 1 
states, which are the only states populated in muon capture on T = targets like 12 C, 16 0, 
and 40 Ca, residing at lower excitation energies. Consequently, the energy of the neutrinos 
E u , which are released after the capture to these states, are slightly higher. As the muon 
capture rate is proportional to E^ 1+x ' (A is the order of the spherical Bessel function of the 
corresponding operator) jjj], it is increased when the residual interaction is weakened. The 
same sequence of arguments shows that the rate is lowered if the interaction is increased. 
These expectations are confirmed by our calculations. For example, the total muon capture 
rate on 40 Ca is changed to 22.77 x 10 5 s _1 (24.18 x 10 5 s _1 ) if the overall strength of the 
Landau-Migdal is increased (decreased) by 10%. When compared to our results, given in 
Table I, we conclude that a 10% variation in the interaction results in a change of the 
total muon capture rate of less than 3%. The same sensitivity is observed for the nuclei 



12 C and 16 0, where the same variation in the interaction changes the rate by less than 3% 



and 2%, respectively. We conclude that the changes in the muon capture rates induced by 
the in-medium renormalization of the form factors qa and gp are noticeably larger than its 
sensitivity to reasonable variations of the residual interaction. 

Before drawing conclusions from our calculations, it is perhaps worthwhile to briefly 
review other calculations of the muon capture rate for the nuclei we are considering. In 



the classical paper of Foldy and Walecka 15], the authors relate the dipole capture rate to 



the experimental photo-absorption cross section. In addition, they use symmetry arguments 
to relate the vector and axial vector nuclear matrix elements. In its slightly more modern 
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version |T[, which uses the present free nucleon coupling constants, the calculation gives total 
capture rates quite close to ours in Table I for 12 C and 16 0, and perhaps 20% higher than 
our numbers (and the experiment) for 40 Ca. 

Another calculation worth mentioning is the treatment of 16 O by Eramzhyan et al. |L6 . 
This calculation employs a truncated shell model with ground state correlations included 
and standard free nucleon coupling constants. For the bound states in 16 N the results are 
similar to ours in Table II, particularly to our entries for the Landau-Migdal force. For the 
transitions to unbound states (which the authors treat as shell model, i.e., bound anyway) 
they obtain a total capture rate of 1.02 x 10 5 s, within 10% of our result. On the other 
hand, Ohtsuka's [17j calculated capture rates in 12 C and 16 are larger than ours and the 
experiment by a substantial factor of 1.5 and 2.5, respectively. The author attributes this 
overestimate to the neglect of ground state correlations in the adopted nuclear model. In 
fact, the same calculation overestimates the photo-reaction cross section, which has nothing 
to do with weak interactions or axial vector current, by a similar factor. 

All of these calculations are therefore basically compatible with our result, and suggest 
that the experimental rate is best reproduced with the free nucleon coupling constants. 
A similar conclusion follows from the comparison of the measured l2 C(i> e , e~) 12 B* cross 
sections for the v e from stopped muon decay ||18|| , with the results calculated within the 
same approach as employed here |19|| . 



In contrast, the recent calculation in Ref. |2(J is performed quite differently. It uses 
an approach based on the local density approximation to the infinite nuclear matter. In 
Ref. p0| , the capture rate is reduced by a factor of about two by the strong nuclear renor- 



malization. This renormalization seems to include both the effects of residual interaction, 
which reduce the rate as we argued above, and the effect of in-medium renormalization. It 
is difficult to separate the two, and therefore difficult to compare our results with those of 
Ref. 0. 

Our calculations show that the continuum RPA method with a standard, unadjusted 
residual interaction describes the muon capture rates in the T = nuclei quite well, pro- 



vided that free nucleon form factors are used. In particular, quenching of the axial current 
coupling constant suggested in the analysis of the Gamow- Teller strength would result in a 
noticeable disagreement with the data. We checked that this conclusion is not sensitive to 
reasonable variation in the residual force. We stress that the muon capture is dominated by 
the transitions to the negative parity dipole and spin-dipole collective states. Our conclusion 
is therefore only relevant for such states. 
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FIGURES 

FIG. 1. Muon capture in 12 C as a function of the neutron energy E(n). Full line is the total 
rate, dashed line is for the 1~ multipole, short dashed line is for the 2~ multipole, and the dotted 
line is for the 2 + multipole. The calculation is for the Landau-Migdal interaction. 

FIG. 2. Muon capture in 16 0. Full line is the capture rate for the Bonn potential and the 
dashed line is for the Landau-Migdal force. 

FIG. 3. Muon capture rate in 40 Ca as a function of the excitation energy in the final nucleus 
40 K. The calculation is for the Landau-Migdal force. 
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TABLES 

TABLE I. Comparison of calculated total muon capture rates with experimental data pj}| . For 
12 C and 16 the capture rates to particle-bound states have been subtracted. The rates are given 

in 10 5 s" 1 . 

Target nucleus Experiment Bonn potential Landau-Migdal potential 

12 C 0.320±0.01 0.342 0.334 

16 0.924±0.01 0.969 0.919 

40 Ca 25.57 ±0.14 26.2 23.4 
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TABLE II. Partial muon capture rates for the bound states in 12 C(/u ,t0 12 B* (upper part) 



and 0(fi~ jU^) N*(lower part) in units of s , calculated for the Landau-Migdal force (LM) 
and the Bonn potential (BP). In each part the theoretical results are shown in the top two lines, 
followed by the measured data. Note, that the values marked with # were assumed to be 0, as the 
+ — > 2 + transition is second-forbidden. 



Source 



u,(l+) 



0/(1") 



0/(2") 



,(2+) 



(LM) 

(BP) 
Ref. fl22l 



Ref. |3| 



Ref. |4| 



Ref. H| 



Ref. |5| 



25400 
22780 
6290±300 

6000±400 
5700±800 
6280±290 



220 

745 
720±175 
1080±125 
890±100 
700±400 
380±100 



40 
25 
10±230 
60±200 
170±240 
400±600 
120±80 



< 1 

< 1 
0# 
0# 
0# 

200±400 
270±100 





Source 


w(O-) 


0/(1-) 


o;(2-) 



(LM) 

(BP) 
Ref. || 
Ref. |7| 
Ref. |8| 
Ref. HI 



1.45 
1.85 
1.1±0.2 
1.6±0.2 
0.85±0.06 
1.56±0.17 



1.75 
3.10 
1.9±0.1 
1.4±0.2 
1.85±0.17 
1.31±0.1 



8.10 
8.65 
6.3±0.7 



8.2±1.2 
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TABLE III. Same as Table I, but calculated for the three renormalization models of gA and 
gp, defined in Eqs. (3-5). The calculations have been performed for the Landau-Migdal force. The 
rates are given in 10 5 s _1 . 

Target nucleus Experiment Model 1 Model 2 Model 3 

12 C 0.320±0.01 0.245 0.252 0.267 

16 0.924±0.01 0.682 0.698 0.736 

40 Ca 25.57±0.14 17.3 17.8 18.8 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9405030vl 



differential capture rate in 1000/(s MeV) 




This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9405030vl 



differential capture rate in 1000/(s MeV) 




This figure "figl-3.png" is available in "png" format from: 
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